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Abstract: 

We investigate the dynamics of the geometric transitions associated to compactified 
spacetimes. By including the dynamics of gravity we are able to follow the evolution 
of collapsing cycles as they attempt to undergo a topology changing transition. Rather 
than achieving this singular geometry we find that one of two scenarios occur, depending 
on the initial conditions. Either a horizon forms, shielding a curvature singularity, or 
the cycle re-expands after an initial contraction phase. For the case where a horizon 
forms we identify the final state with a known analytic black-hole solution. We also 
show use our results to demonstate a novel compactification mechanism, owing to the 
asymptotic structure of this black-hole solution. 
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1. Introduction 

The term string-theory can be interpreted in a number of different ways, including 
26-dimensional bosonic string theory, 11-dimensional M-theory and 10- dimensional su- 
perstring theory. The common element is that these theories exist in dimensions other 
four. Given that we only experience four dimensions we are forced, if string-theory 
is correct, to explain why we do not observe these extra dimensions. There are two 
known mechanisms to render extra dimensions unobservable, compactification U and 
brane- worlds 0; we shall be studying compactification, whereby the full spacetime is 
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constructed out of four large dimensions and some small compact manifold. If we use 
the framework of superstring theory then the compact manifold Ai has six dimensions, 
and if we require the low-energy theory to be supersymmetric this imposes certain re- 
strictions on the type of manifold Ai can be. In this paper we shall not be switching on 
any of the fluxes sourced by D-branes and therefore we find that Ai must have special 
holonomy ||, implying that Ai is a Calabi-Yau space. A problem now arises in the 
lack of a unique way to choose Ai, there exist many different Calabi-Yau manifolds 
on which the compactification can take place. Some only differ by parameters, giving 
a continuous degeneracy, but some have completely different topologies, leading to a 
discrete degeneracy. Calabi-Yau manifolds with the same topology are grouped into 
moduli spaces, with the spectrum of massless moduli fields in the low-energy theory 
depending solely on the topology of Ai. 

This situation of having to choose one Calabi-Yau from the many was made more 
palatable when it was realised that many of these Calabi-Yau manifolds were in fact 
connected by finite length paths in moduli space || |7|, |8[], with the connecting geometry 
being singular. One way to picture these transitions between manifolds is to study the 
cycles within them, for example it may be that certain cycles collapse to zero size on 
one side of the transition and expand as different cycles on the other. During such 
a transition the Hodge numbers need not change, giving only a change of intersection 
numbers, such as a flop transition where an S 2 collapses and re-expands as a different S 2 . 
Even though the Hodge numbers are unaffected this still constitutes a topology change 
and so is necessarily singular from the geometrical perspective. Remarkably, string 
theory is able to make sense of these singular geometries by the appearance of new 
light states corresponding to D-branes wrapping the collapsing cycles; the dynamics 
of the low-energy theory of flop transitions has been studied in ||, [H], |TT[. Other 
transitions may change the Hodge numbers, such as a conifold transition [f|, || where an 
S 2 collapses and reappears as an S 3 , see also |J. Again, such a singular transformation 
can be made regular within string theory ||12|| , and the corresponding low-energy theory 



may be studied |H| along with its dynamics |14], [L5| . 

We have just described one of the key uses of singular geometries within string- 
theory, namely they connect together different moduli spaces. This was achieved by 
the appearance of extra light states coming from the wrapped D-branes around the 
collapsed cycles. These "extra" states also play a crucial role in building low-energy 
models with chiral fermions and non-Abelian gauge fields, allowing one to evade a no-go 



theorem [16| for low-energy chiral fermions coming from compactification [17, 18, 19, 



The types of singularity that have proven useful in the above constructions are 
conical, taking the local description of a discrete quotient of a smooth manifold X, 
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X = X/T, where T is a finite symmetry group; the singularity is then the fixed point 
set of r. Of particular interest are those singularities which allow a smooth resolution, 
by blowing up certain cycles of zero size contained within the fixed point set of T. This 
is achieved in practise by a surgery which replaces a ball around the conical singularity 
with a ball of a smooth special-holonomy space [[30], |31], E3|| . The moduli of this special- 
holonomy space then appear in the low-energy theory as moduli fields [E3]. 



In this paper we study the dynamics of the resolved spaces while the cycles are 
collapsing. Unlike previous studies on the dynamics of such spaces ||, [H| |y], [14], [15 



we are particularly interested in the gravitational properties of collapsing cycles; more 
specifically, whether a horizon forms as the cycle becomes small. Should a horizon 
form in the higher dimensional theory, then this would render the low-energy theory 
based on the moduli fields inapplicable near the conical singularity, implying that the 
dynamics of topology changing processes is more complicated than simply studying the 
dynamics of the moduli fields in the low-energy description. 

To have a specific model we shall be studying the resolution of the C 2 /Z2 orbifold. 
Singularities of this type are common in the construction of compact manifolds of 
special holonomy f24|, where one typically starts with a torus, imposes some reflection 



symmetries, and then replaces the fixed points of the symmetries smooth manifolds. In 



our case the singularity is resolved with the Eguch-Hanson instanton pi] , |22| . We begin 
our study in section |2] with a discussion of four- dimensional Euclidean instantons, and 
introduce the Eguchi-Hanson solution along with its size modulus for the collapsing 
S 2 . In section || we describe the method we use to simulate Einstein's equations, with 
the initial data presented in section |j. We finish by giving our results in section [5] and 
conclusions at the end. 

2. Gravitational instantons 

To have an explicit construction of a compact manifold with special holonomy using 
the methods of Joyce |24j] one needs to have explicit Ricci-flat metrics of special holon- 



omy; a number of such metrics are known for non-compact manifolds [EI, 26, E_7], 28 
One then replaces a neighbourhood of the singularity with a region of the non-compact 
space, using some form of smoothing at the join, to construct the metric on the smooth 
compact manifold. As we are interested only in the dynamics near the putative con- 
ical singularity we only sudy the dynamics of the non-compact space, assuming that 
information far from the region of interest has no impact over the timescales of our 
simulations. 

The singularity we are interested in resolving is the simplest case, C 2 /Z2, and for 
that we need to know something about Ricci-flat metrics on four-dimensional Rieman- 
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nian manifolds, also known as gravitational instantons. 



2.1 Eguchi Hanson metric 

The Eguchi-Hanson metric is a regular self-dual, hyper-Kahler metric in four-dimensions 
and has the asymptotic structure of C 2 /Z 2 |2l|, |22|, i.e. it is a resolution of the C 2 /Z 2 
conical singularity. It is constructed as a cohomogeneity-one metric with squashed 
three-spheres as the level surfaces and has the explicit form 

ds 2 EH (l) = a(p)~Mp 2 + [(af + al) + a(p)a 2 ] , (2.1) 

a(p) = 1 -(-) . (2.2) 



We have used the conventional left-invariant one-forms of SU(2) which satisfy 

da; = -\eijkOj A cr fc , (2.3) 

and the parameter I is a constant parameter i.e. a modulus of the solution. 

From the above form of the metric we see that there is an apparent singularity at 
p = I, we get a clearer understanding of its nature if we look close to this region using 
the following co-ordinates, 

p = l + T- (2 - 4) 

This results in the metric taking the form 

ds 2 = [d# + #o*]+£tf + o*), (2.5) 

which clearly shows that the apperent singularity at R=0 (p = I) is just a co-ordinate 
artefact, and that the manifold looks locally like a product of flat space and a two- 
sphere of radius 1/2; this type of removable singularity is termed a bolt singularity p5| . 
It is the finite size of this two-sphere which has resolved the singularity, by taking / 
to zero in (|2.1|) we explictly see the metric becomes C 2 /Z 2 . (The Z 2 comes from an 



identification required to make the origin of the resolved space regular ||22|| .) 
2.2 Moduli evolution. 

In our simulations we shall be working in five dimensions, with the four spatial dimen- 
sions initially taking the form of the Eguchi-Hanson metric. This is an exact solution 
of Einstein's equations, so to get the spacetime to evolve we must give the metric some 
momentum. Before presenting the full numerical approach we should see what we can 
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find analytically. As pointed out by Manton in the context of BPS monopoles |29j one 
can understand the low-energy dynamics of a system by allowing its moduli to have 
a small time dependence. We introduce a time-dependent modulus L(t), such that 
L{t = 0) = I. In our case we can use the Einstein-Hilbert action 



S = J dt d 4 x y/=gR, (2.6) 

with the metric 

ds 2 = -dt 2 + ds 2 EH (L(t)) , (2.7) 
to derive the effective action for the modulus L(t), 



n 2 f . f d 



= t J dt [at 12 ) • ( 2 - 8 ) 

This leads to the conclusion that the moduli space approximation predicts that L(t) 2 
evolves linearly with time. Also, given that the boundary of the moduli space is L = 
we see that the modulus L can reach the boundary in finite time. 

3. Numerical Simulation 

The idea now is to start with a regular five-dimensional metric and evolve it numerically 
using Einstein's equations. A similar project was undertaken by Bizon el al |32| who 



realised that Birkhoff's theorem could be evaded in five spacetime dimensions. Their 
initial metric however contained a nut singularity at the origin as opposed to the bolt 
in ([2.7]) . While both nut and bolt singularities are removable co-ordinate singularities 



and so are regular analytically, we found that the numerical techniques of [|32| broke 



down at the origin. After numerous attempts, using different forms for the metric, we 



settled on the ADM formalism for evolving the metric |33|, |34]]. This involves three 
metric functions along with their momenta, 

ds 2 = -dt 2 + a{t, r) 2 dr 2 + b(t, r) 2 (a\ + o\) + c{t, r)V 2 . (3.1) 

The ADM equations then give first order evolution equations for the metric functions 
and their momenta, along with a Hamiltonian constraint and momentum constraints. 
Note that we have chosen the gauge such that the lapse is unity and shift vector 
vanishes, we have also enforced the bi-axial form of the Eguchi-Hanson metric for 
simplicity. It is certainly possible for ( ^.1| ) to describe the metric we are interested 
in, ( p.7|) , however if we were to simply equate the r co-ordinate of ( p.l|) with the p 
co-ordinate of (|2.7|) we would find that a(t, r) diverges at the origin. This clearly 
causes unacceptable numerical problems, so we need to choose our parametrization 
more carefully, as desribed in the next section. 
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3.1 Metric 

To aid the stability of the algorithm, particularly in establishing sensible boundary 
conditions, we require that 

1. All three variables and all three momenta to remain even and finite at the origin. 

2. All three variables and all three momenta to tend to finite (maybe zero) values 
asymptotically. 

3. We minimize the amount of division by variables in all equations of motion. 



To that end we evolved the following form for the metric 



ds 2 = -dt 2 + 



v \ 2 

1 "T 



e 1A Ar' + 1' 



r \ 4 

1 '7 



<Y* \ 2 

1 <7 



e 2C o\ 
(3.2) 

In order to impose all the necessary boundary conditions at once, keeping the equa- 
tions regular at the origin, we used techniques outlined in |39], [40| . This involved the 
introduction of three new variables Da,D b and D c , to replace the spatial derivatives, 



D A 
D c 



B' + 
C' + 



4r 



2r 3 
r 



P + r 2 

We also introduced the momenta Ka,K b and K c , defined as 



(3.3) 

(3.4) 
(3.5) 



K A 
K b 
K c 



-A, 
-B, 
-C, 



(3.6) 
(3.7) 
(3.8) 



where ' indicates derivative with respect to time and ' indicates derivative with respect 
to r. The Di were chosen so as to be odd at the origin as these were simple boundary 
conditions to impose. Actually, the full set of boundary conditions at the origin may 
be found by requiring local flatness 1311 , in which case we find 



A(t,r) ~ A (t) + O(r 2 ), 
D A {t,r) ~ 0(r), 
K A (t,r)~K A (t)+O(r 2 ), 



(3.9) 
(3.10) 
(3.11) 
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with similar relations for the functions B(t,r), C(t,r), Ks(t,r), Kc(t,r) at the origin. 
We also find that 

A°(t) = C°(t), K° A = K° C . (3.12) 

By giving the metric functions some initial momentum the spatial part of the metric will 
cease to remain Eguchi-Hanson, however it will retain some Eguchi-Hanson features, 
at least for early times. Notably, the bolt singularity at the origin will remain, still 
describing a two-sphere of radius L(t)/2. We used the value of B at the origin to define 
this L(t) at later times. 

L(t) = 2exp(B°(t)) (3.13) 
Note that for L(t) to vanish, then B(t,r = 0) must diverge. 



4. Initial conditions 



In the parametrization of (3.2p, using the co-ordinate r — R of ( |2.4j ) (|2.1|) we find that 
our initial conditions for the metric functions take the form 



A = - In 
2 



"i 2 



B 



2 
1 



In 



1 - 

r r \ 2 
1 



-i 



I 2 + r 2 



j u 7 



1 + 



-i 



C = -ln 



4r 2 




(4.1) 



If we impose vanishing momenta then this would constitute an exact solution of the 
equations of motion. As a check of our numerics we do indeed find that the system 
remains static. As we want to evolve the Eguchi-Hanson metric toward the conical 
singularity we must impose some non-vanishing momentum for the metric functions. 
This is not a completely trivial task given that general relativity imposes constraints 
coming from the gauge fixing (appendix |A]). The two constraints, Hamiltonian and 
momentum, mean that once A(t = 0, r), B(t = 0, r) and C(t — 0, r) are fixed according 
to ( ft-lp there is one free function left to describe the momentum. To fix this function 



we take our motivation from the moduli space approximation of section [2.2| and find 
that initially we have 

K B = -B = -k (t^l ), (4-2) 



P + r 2 
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so we are able to choose an L and derive from this Kb- we imposed that L was required 
to be: 



1. even at the origin. 

2. finite and negative at the origin. 

3. exactly zero far from the origin. 

4. continuous and differentiable to first order. 

The first condition ensures that Kb is even, the second means that we push the Eguchi- 
Hanson space towards the conical singularity. The third condition is imposed so that 
only the form near the origin is important, and that the non-compact nature of Eguch- 
Hanson does not affect the evolution. The final condition gives a smooth profile for us 
to evolve. 

Only one of the three momenta, Kb-, was specified explicitly by L, with the other 
two being derived from the constraints (|A.l , |A2| ) using a 4th order Runge-Kutta algo- 
rithm. 

We decided on an L, taking the form: 

-L (l~& 2 r<rO 

r > rO K ' J 

where Lq is a positive constant (the magnitude of L at the origin) and rO is another 
constant which determines the outer radius of the non-zero L. With the initial data 
fixed we now evolve the system according to the equations of motion laid out in ap- 
pendix [TJ. To acheive this we used a 4th order Runge-Kutta algorithm, while keeping 
a check that the constraints of appendix [A] remained small; typically they were of order 
0.005. 



5. Results 

5.1 Apparent horizons 

The outgoing radial null geodesies can usually be seen to be increasing in area, but 
after the formation of the apparent horizon, which is a null trapped surface, the null 
geodesies are no longer increasing in area all null rays are in fact converging. This 
apparent horizon shows that there exist points in space for which all null geodesies are 
unable to diverge to X + . These points are therefore behind an event horizon. This 
event horizon must include everything within the apparent horizon (and maybe more), 
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however the apparent horizon is far easier to detect and it proves the existence of 



the event horizon ||35|| . In order to determine whether a black hole has formed we 
continually looked for such an apparent horizon. 

At the apparent horizon, a congruence of null geodesies no longer increase their 
area and we may write, 







dArea 

dt 
d Area 

dt 



null 
+ 



d Area 




dr 


dr 


t 


At. 



null 



The area is determined by Area ~ b 2 c in terms of Q), with g ] 
also be calculated in terms of the extrinsic curvatures [BH| to give 







-2K B -K C + 



2D B + D a + h 



l + 4(i) V 



(5.1) 
(5.2) 

1/a. This may 
(5.3) 



at an apparent horizon. 

5.2 Formation of black holes 

Depending upon our input parameters, L and rO, there were three possible outcomes 
to adding the momentum: 

1. For sufficiently low L and rO, there was insufficient initial momentum to observe 
the creation of either a black hole or a singular topology. 

2. For an intermediate range in the parameters the system produced an apparent 
horizon. After initially increasing, the area of the apparent horizon converged to 
a constant value. 

3. For large initial parameters the system already contained an apparent horizon 
simply due to the initial conditions. 

Our interest is directed to case 2. In this case a black hole forms, which was not 
initially present. 

Before we present the results for a range of parameters we focus on a single example 
where we took rO = 1.0, L Q = 2.3. In Fig. |l]a we plot the time dependence of L as 
measured at the origin, and also give the area of the horizon - both in units of /. What 
the figure shows is that L is monotonically decreasing, with L 2 decreasing approximately 
linearly in the initial phase. This is consistent with the expectations from the moduli 
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space approximation of section [2.2| . However, at some point an apparent horizon forms 
(t ~ 0.11 in the simulation) which then implies an event horizon exists - something the 
moduli space approximation does not account for. At this point we can no longer trust 
any low-energy dynamics derived from the moduli space approximation. We also see 
from Fig. |TJa that the area of the apparent horizon increases initially, but the settles 
down to a fixed value. Presumably this corresponds to the formation of what would 



become a static black hole; we shall discuss this further in section |5.3j The figure also 
shows that L, as defined at the origin, reaches zero in finite co-ordinate time t (t ~ 0.3 
in the simulation). This corresponds to a divergence in the metric function B(t,r) and 
is in fact a curvature singularity. Fortunately this is hidden behind the horizon. 

In order to get a clearer understanding of the causal structure of our solution we 
present in Fig. [TJh> a plot showing various radial outgoing null geodesies. Superimposed 
on this is the curve showing the location of the apparent horizon. We see that initially 
the null rays continue outwards and, given the asymptotically locally flat structure of 
Eguchi-Hanson, reach null infinity. However, some time later the outgoing null rays 
near the origin turn around and head towards r = 0. The presence of such null rays 
indicates that a horizon has formed, and this is confirmed by the existence of the 
apparent horizon. 
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Figure 1: (a) The steady decrease of I? and the formation of the apparent horizon, (b) The 
outgoing radial null geodesies and the apparent horizon which later formed. 
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5.3 The resulting black hole 

By the time the program ends (due to the curvature singularity at r = 0) the apparent 
horizon has settled to a single area which can be measured. The natural question 
is "what is the final state?" Given that we cannot run the simulations beyond the 
curvature singularity we can only offer a conjecture to answer this question. However, 
given that the horizon has converged to a constant value we believe that it is reasonable 
to suggest a five-dimensional black-hole is in the process of forming. The black-hole 
which fits our requirements was written down in its Kaluza-Klein dimensionally reduced 
form in |37, |38|. Written in its five- dimensional form this black hole looks like |36[ 



ds 2 = -fdt" + jdr* + £ [k{a\ + al) + a 3 2 ] , (5.4) 



f ^,2 



fir) = y 3 + (5-5) 



( ~,2 ^f»-*2 

k (r) = \Z zkr (5-6) 



I /T"2 /y>2 |2 

V co 



and describes a static black-hole with a squashed three-sphere for a horizon at r = r + . 
The radial co-ordinate range is < r < and the parameter range is < r + < r m . 
If we accept that this is the end state of the Eguchi-Hanson collapse then we are free 
to evaluate the squashing function k(r) at the horizon, provided it too has settled to a 
single value before the program's end. 

The asymptotic structure of the black-hole is interesting in that it is not asymp- 
totically flat, rather it is asymptotically locally flat and takes the form [ ]5E|] 

2 

ds 2 = -dT 2 + di? 2 + R 2 dQ S 2 + r -^ x 2 - (5-7) 

So, locally this looks like MS 1 ' 3 ' x S 1 , where the circle has radius we can find 

the parameters r + and by evaluating the area of the horizon, and the squashing 
parameter on the horizon ( k (r + ) = k + ), 

(5.8) 
(5.9) 

This result gives us a rather novel method for dynamical compactification. Suppose 
that instead of starting with a compact manifold, where a portion of Eguchi-Hanson 
space has been glued in, we start with the full Eguchi-Hanson space with its four "large" 
spatial dimensions. Then our results show that this evolves to a space where one of 
the spatial dimensions compactifies to a circle, giving three "large" dimensions and one 
"small". 
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5.4 Varying the initial data. 



The black hole's area and the extent to which the angular part is squashed depends 
on the initial conditions; in our parametrization ( |4.3| ) this means changing L and rO. 
The results of varying L while keeping rO constant are shown in Fig. 0. 

As described earlier, for values of L which are too small the system never produces 
any sort of horizon, no singularity is formed and the value of L drops for a small 
amount then begins to rise again, there is not enough energy to form a black hole or a 
singularity. 

Alternatively, if we take L to be too large then the initial data already contains an 
apparent horizon, rather than forming one dynamically. The results we present in Fig. 
^| cover the intermediate range where there is enough localized energy to form a black 
hole, but not so much that it is there at the start of the simulation. 

In Fig. 0, where rO = 1.0, the apparent horizon forms dynamically for 0.7 < Lq < 
2.8 and it's area can be seen to converge and be measured. Over the duration of the 
simulation the squashing parameter was seen to converge for the range 2.0 < L < 2.8, 
and in all these cases it converges to a value greater than one. This is consistent with 
the numerical squashing parameter being identified with the analytic form of k + , given 
in (|5.6|), which must also remain greater than one at the horizon. 



1 5 1 ■ ■ ■ ■ ■ ■ 1 2 




Figure 2: The effect on the final area and extent of squashing (only for values for which it 
has converged) at the horizon due to differing the initial Lq (rO = 1.0). 



If we vary both the values of Lq and rO we can find the resulting area in a great 
many cases, this is shown in Fig. [| The range in which a horizon forms at a late time 
is shown and given a shading scale to indicate the area of that horizon. Within the 
region marked A, there is insufficient energy to form a horizon at all. Region B marks 
the existence of a horizon within the initial conditions. 



Area 




Figure 3: The effects on the creation of a black hole and it's final area due to differing the 
initial Lq and rO. 

6. Conclusions 

The possibility of transitions which change the topology of a Calabi Yau manifold, 
is intriguing. It suggests a method of connecting, seemingly discrete, moduli spaces 
of Calabi Yau manifolds by their common singularities. If the transitions can occur 
dynamically, then the topology of the compactified dimensions of string theory may 
change in time. Based on this expectation one can study the low-energy dynamics 
coming from such a transition. 

To study these topology changing processes we perform numerical simulations of 
an Eguchi-Hanson spacetime with a collapsing two-cycle. Our results highlight the 
importance of gravity during these events where cycles are collapsing, showing that 
either horizons form during the process, or the cycle re-expands. In either case one 
concludes that the gravitational effects prevent the naive collapse of the cycle. 
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We have presented evidenvce that, in the case where a horizon forms, the final state 
of the evolution is a black hole, where the horizon is a squashed three-sphere |37], |38 



Such a black-hole has an interesting asymptotic structure, namely there is a compact 
circle at infinity, and this leads us to an unexpected mechanism for compactificaction. 
If, instead of picturing the Eguchi-Hanson space as a portion of a compact internal 
space, we start with the full Eguchi-Hanson space, with its four "large" dimensions, we 
see that the final state has a compact dimension and corresponds to the Kaluza-Klein 



black- hole of [37 
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Appendices 
A. Constraints 

The metric produced the following constraints, equations the variables must always 
conform to. These were imposed as initial conditions and later monitored to test the 
program's accuracy. 

The Hamiltonian constraint: 

0= -K A (2K B + K c )-K B (K B + 2Ko) + ^-^ (A.l) 

+ \ ( -D A (D c + -) + 2 D B (D c + -) - 2 D A D B + 3 D% + 2 D' B + D' c + ^ + D 2 C ) 
a z \ r r r J 

Also the momentum constraint: 

K A {D c + 2D B + l -)=2K' B + K' c + 2K B D B + K c {D c + ^-). (A.2) 

where: 

a 2 = ( 1+4 (9 2 ) b 2 = 4/ 2 (l + (y) 4 ) c 2 = 4r 2 (l + 

(A.3) 

Apparently singular terms within these constraints did not produce any instabilities 
as they do not feed back into the equations used to evolve the system, they were only 
used for testing purposes. 



- 14 - 



B. Equations of motion 

We found the equations of motion from the field equations using the ADM formalism 



33| [34] . Also we added multiples of the momentum and Hamiltonian constraints to 
remove as many potentially singular terms from our equations of motion. This resulted 
in the following equations of motion. 

A = -K A 
B = -K B 
C = -K c 



K A = + 2D B (D c + -)) + ^ - -J + K\ - K% - 2K B K C 



K B = \{D B D a -2D 2 b -D' b -Db{Dc + -))-^ + ^ + 2K 2 b + KbK c + K b K a 
a z ' rob 

K c = ±(3D%-2D b D a + 2D' b ) + ^-^-K b + K 2 c -2K b K a 



D A = 


~K' A 


D B = 


-K'b 


D c = 


-K'c 



fB.r 



Where a, b and c are defined in ( 14 .3| ). The only potentially singular term remaining 



(which could have produced instabilities) is D B /r, analytically this is regular as D B is 
odd. Numerically it was sufficiently stable to allow the program to run its course. 

References 

[1] T. Kaluza, Sitzungsber. Preuss. Akad. Wiss. Berlin (Math. Phys. ) 1921, 966 (1921). 

[2] V. A. Rubakov and M. E. Shaposhnikov, Phys. Lett. B 125, 136 (1983). 

[3] P. Candelas and D. J. Raine, Nucl. Phys. B 248, 415 (1984). 

[4] P. Candelas and X. C. de la Ossa, Nucl. Phys. B 342, 246 (1990). 

[5] R. Gwyn and A. Knauf, arXiv:hep-th/0703289. 

[6] P. Candelas, P. S. Green and T. Hubsch, Phys. Rev. Lett. 62, 1956 (1989). 

[7] P. Candelas, P. S. Green and T. Hubsch, Nucl. Phys. B 330, 49 (1990). 



- 15 - 



[8] B. R. Greene, arXiv:hep-th/9702155. 

[9] M. Brandle and A. Lukas, Phys. Rev. D 68, 024030 (2003) [arXiv:hep-th/0212263]. 

[10] L. Jarv, T. Mohaupt and F. Saueressig, JHEP 0312, 047 (2003) 
[arXiv:hep-th/0310173]. 

[11] L. Jarv, T. Mohaupt and F. Saueressig, JCAP 0402, 012 (2004) 
[arXiv:hep-th/0310174]. 

[12] A. Strominger, Nucl. Phys. B 451, 96 (1995) [arXiv:hep-th/9504090]. 

[13] B. R. Greene, D. R. Morrison and C. Vafa, Nucl. Phys. B 481, 513 (1996) 
[arXiv:hep-th /9608039] . 

[14] A. Lukas, E. Palti and P. M. Saffin, Phys. Rev. D 71, 066001 (2005) 
[arXiv:hep-th/0411033]. 

[15] E. Palti, P. Saffin and J. Urrestilla, JHEP 0603, 029 (2006) [arXiv:hep-th/05 10269]. 

[16] E. Witten, "Fermion Quantum Numbers In Kaluza-Klein Theory," reprinted in 
A. Salam and E. Sezgin, Supergravities in diverse dimensions Vol 2, Amsterdam, 
Netherlands: North-Holland (1989). 

[17] B. S. Acharya, arXiv:hep-th/0011089. 

[18] M. Atiyah, J. M. Maldacena and C. Vafa, J. Math. Phys. 42, 3209 (2001) 
[arXiv:hep-th/0011256]. 

[19] M. Atiyah and E. Witten, Adv. Theor. Math. Phys. 6, 1 (2003) [arXiv:hep-th/0107177]. 

[20] B. Acharya and E. Witten, arXiv:hep-th/0109152. 

[21] T. Eguchi and A. J. Hanson, Phys. Lett. B 74, 249 (1978). 

[22] T. Eguchi, P. B. Gilkey and A. J. Hanson, Phys. Rept. 66, 213 (1980). 

[23] A. Lukas and S. Morris, Phys. Rev. D 69, 066003 (2004) [arXiv:hep-th/0305078]. 

[24] D. Joyce, "Compact Manifolds with Special Holonomy," Oxford University Press, 2000 

[25] G. W. Gibbons and S. W. Hawking, Commun. Math. Phys. 66, 291 (1979). 

[26] G. W. Gibbons and S. W. Hawking, Phys. Lett. B 78, 430 (1978). 

[27] G. W. Gibbons, D. N. Page and C. N. Pope, Commun. Math. Phys. 127, 529 (1990). 

[28] M. Cvetic, G. W. Gibbons, H. Lu and C. N. Pope, Phys. Rev. D 65, 106004 (2002) 
[arXiv:hep-th/0108245]. 



- 16 - 



[29] N. S. Manton, Phys. Lett. B 110, 54 (1982). 

[30] G. W. Gibbons and C. N. Pope, Commun. Math. Phys. 66, 267 (1979). 
[31] D. N. Page, Phys. Lett. B 80, 55 (1978). 

[32] P. Bizon, T. Chmaj and B. G. Schmidt, Phys. Rev. Lett. 95, 071102 (2005) 
[arXiv:gr-qc/0506074]. 

[33] R. Arnowitt, S. Descr and C. W. Misner, arXiv:gr-qc/0405109. 

[34] B. Kelly, P. Laguna, K. Lockitch, J. Pullin, E. Schnetter, D. Shoemaker and M. Tiglio, 
Phys. Rev. D 64, 084013 (2001) [arXiv:gr-qc/0103099]. 

[35] J. Thornburg, Phys. Rev. D 54, 4899 (1996) [arXiv:gr-qc/9508014]. 

[36] H. Ishihara and K. Matsuno, Prog. Theor. Phys. 116, 417 (2006) 
[arXiv:hep-th/0510094]. 

[37] G. W. Gibbons and D. L. Wiltshire, Annals Phys. 167, 201 (1986) [Erratum-ibid. 176, 
393 (1987)]. 

[38] P. Dobiasch and D. Maison, Gen. Rel. Grav. 14, 231 (1982). 

[39] M. Alcubierre and J. A. Gonzalez, Comput. Phys. Commun. 167, 76 (2005) 
[arXiv:gr-qc/0401113]. 

[40] M. Ruiz, M. Alcubierre and D. Nunez, arXiv:0706.0923 [gr-qc]. 



-17- 



